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The LHCb Collaboration at the Large Hadron Collider at CERN discovered two pentaquark states 
and Rc( 4450). These two hidden-charm states are interpreted as the loosely bound Ec(2455)Z)* and E*(2520)D* 
molecular states in the boson exchange interaction model, which provides an explanation for why the experimen¬ 
tal width of R(,(4450) is much narrower than that of Rc(4380). The discovery of the new resonances Rc(4380) 
and Rc(4450), indeed, opens a new page for hadron physics. The partners of Rc(4380) and Rc(4450) should be 
pursued in future experiments. 

PACS numbers: 14.20.Pt, 12.39.Pn 


In the pioneer paper [1], Gell-Mann indicated that the mul¬ 
tiquark states should exist along with the simplest structures 
for baryons which are composed of three valence quarks and 
mesons which contain a quark and an antiquark. Indeed, from 
the point of either mathematics or physics, there is nothing to 
forbid the existence of such exotic states. In quantum mechan¬ 
ics, the multiquark states are nothing new but just higher Fock 
states. In gauge field theory, the QCD principle allows the 
existence of multiquark states and hybrids which contain not 
only quarks but also gluonic degrees of freedom. The mul¬ 
tiquark states might be suppressed by higher orders of as in 
perturbative QCD. However, the confinement is a totally non- 
perturbative QCD effects. Such suppression is not so drastic. 
Unless nature prefers the simplest structures, one has reason 
to believe that the multiquark states should exist. 

From the aspects of both theory and experiment, the ground 
state mesons and baryons are well arranged in the simplest 
structures; singlets and octets for mesons, octets and decuplets 
for baryons. No flavor exotic states have been observed or 
required to exist. However, for the mesons whose masses are 
above 1 GeV, the mixtures among regular quark structures and 
hybrids, glueballs, and multiquark states should be non-zero 
in order to fit the spectra [2, 3]. 

In 2003, the LEPS Collaboration reported the evidence of 
the strange pentaquark state with the quark content uudds and 
very narrow width [4]. Unfortunately this flavor exotic state 
was not confirmed by subsequent experiments [5, 6]. In fact, 
the possible theoretical arguments for the non-existence of the 
stable strange pentaquarks were given in Refs. [7, 8]. Up to 
now, no stable light flavor pentaquark states have been found, 
but light flavor baryons probably have significant pentaquark 
components [9]. 

In the past decade, the researches have rounded the cor¬ 
ner. Many mesonic X, Y, Z particles have been observed at 


Belle, BaBar, BESIII and LHCb. Some of them are identified 
as exotic candidates because they cannot be accommodated in 
the regular qq structures. Corresponding review papers can 
be found in recent literature (for example, see a recent review 
in Ref. [10]). A common point is that all those exotic states 
which are interpreted as molecular states, tetraquarks or even 
four quarks in an anarchy state, contain heavy quarks and an¬ 
tiquarks. 

In other words, the heavy quarks play an important role 
to in stabilizing the multiquark systems, just as in the case 
of the hydrogen molecule in QED. There were the theoreti¬ 
cal predictions of hidden-charm pentaquarks [11-14], espe¬ 
cially the possible hidden-charm molecular baryons with the 
components of an anti-charmed meson and a charmed baryon 
which were investigated systematically within the one boson 
exchange model in Ref. [12]. 

Because of the existence of heavy quarks, the multiquark 
exotic states can survive and light exotic states must mix with 
the regular ones, instead [15]. This conjecture is consistent 
with the fact that all the exotic states which have been ex¬ 
perimentally observed contain hidden charm or bottom. Of 
course, if the allegation that the heavy components stabilize 
multiquark systems, is valid, one suggestion would be natu¬ 
ral, that we should persuade our experimental colleagues to 
search for exotic states which contain open charm or bottom, 
for example, bcqq' etc. Another trend is also obvious, that 
one should be convinced that even though the light pentaquark 
does not individually exist, a pentaquark containing heavy fla¬ 
vors would be possible and has a large probability to be ob¬ 
served in sophisticated experimental facilities. 

Very recently, the LHCb Collaboration observed two res¬ 
onance structures Pc(4380) and Rc(4450) in the T/i/'p invari¬ 
ant mass spectrum of A* —» JlxjipK [16], and their resonance 
parameters are M/)^( 4380 ) = 4380 + 8 + 29MeV, Fp^(4380) = 


2 


205 + 18 + 86MeV, M/>^( 4450 ) = 4449.8 + 1.7 + 2.5 MeV, and 
r/)^( 4450 ) = 39 + 5 + 19 MeV. According to their final state 
y/i/f/7, we conclude that the two observed must not be an 
isosinglet state, and the two Pc states contain hidden-charm 
quantum numbers. A more important feature of these two Pc 
states is that Pc(4380) and 7’c(4450) are close to the thresholds 
of 2c(2455)D* and S*(2520)D*, respectively. 

In this Letter, we propose that the novel resonances re¬ 
ported by the LHCb Collaboration can be identified as the 
exotic hidden-charm pentaquarks with the 2^(2455)/)* and 
I,*(2520)D* molecular configurations. For the Sc(2455)D* 
and Zc(2520)D* systems, the flavor wave functions \I, I^) are 
defined as 

11 1 \ _ [2 |Yt*)++ 1 |y(*)+ n*t)\ 

I 2 ’ 2 / - V 3 ^ ^ I (j) 
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where only these flavors of wave functions with isospin I — 
1/2 match the discussed Pc(4380) and Pc(4450). In the fol¬ 
lowing, we perform a dynamical calculation of the structures 
of 2c(2455)Z)* and S*(2520)Z)* where the constituents inter¬ 
act via one pion exchange (OPE) [12, 26, 27]. The solution 
might help to confirm whether or not bound states for these 
S-wave Sc.(2455)Z)* and S*(2520)Z)* systems exist. To estab¬ 
lish the effective potential which is responsible for binding the 
constituents, we adopt the following effective Lagrangians 


Zp = igTr 




(Q) 




( 3 ) 
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which are constructed under the heavy quark limit and chi¬ 
ral symmetry [17-22]. Here, the notation - 

with v = (1,0) stands for the multiplet field com¬ 
posed of the pseudoscalar P and vector with = 

And the next superfield is composed of 

spinor operators as where 

the notations Sg and Sg are defined as multiplets which, re¬ 
spectively, correspond to - 1/2^ and j’’ - 3/2^ in 6/? 
flavor representations. The axial current is defined as 
A^ = with ^ = exp(/P//;r), and the pion 

decay constant/;r = 132 MeV is taken. Additionally, the ma¬ 
trices P, Sg, and Sg read as 


V2 


V2; 
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Expanding Eqs. (3) and (4), further, we can further get the 
direct effective Lagrangian which will be applied to our later 


calculation, i.e.. 


Zd’D'F - 
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Xs-S-P = 
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where g - 0.59 + 0.07 + 0.01 is extracted from the width of D* 
[23] as is done in Ref. [24], and gi - 0.94 was fixed in Refs. 
[ 12 , 22 ]. 

With the above preparation, the OPE-based potentials for 
the Xc(2455)Z)* and E*(2520)D* systems are deduced, which 
can be related to the scattering amplitude of the l,cD* —> 'ZcD* 
and S*D* —> E*D* processes under adopting the Breit ap¬ 
proximation and performing the Eourier transformation [25], 
where the monopole form factor p'{q~, mp) - (A^ - ot^)/(A^ - 
q^) is introduced to compensate the off shell effect of the ex¬ 
changed meson and describe the structure effect of each effec¬ 
tive vertex. In this form factor, there is a phenomenological 
parameter A which plays an equivalent role as the cut-off in 
the Pauli-Villas renormalization scheme and must be fixed by 
fitting data. We discuss this, later. Einally, the general expres¬ 
sions of effective potentials for the Y.cD* and E*D* systems 
are 

= ^^v2y(A,m,,r)J-o^o, (8) 

^ Jn 

Vi;p.(r) = \^y^Y{K,m„,r)J,gu (9) 

^ Jn 

respectively, where the T(A, m, r) function is defined as 
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In Eqs. (8) and (9), coefficients and Qi {i = 0,1) are related 
to the isospin and quantum numbers of the concerned 

systems. We list them in Table I. 


TABLE I: The values of the J7i and coefficients. Here, S , L, and 
J denote the spin, orbital, and total angular quantum numbers, re¬ 
spectively. S denotes L = 1 since we are interested in the S-wave 
interaction of the Sr(2455)5* and E*(2520)D* systems. 
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By solving the Schrodinger equation with the obtained ef¬ 
fective potentials [28, 29], we can reproduce the masses of 
Pc(4380) and Pc(4450) as shown in Eigs. 1 (a) and 1 (b), 
which supports the allegation that Pc(4380) and Pc(4450) are 
hidden-charm molecular states l,cD* with {I - 112, J - 3/2) 
and Z*Z)* with (I - 112, J - 5/2), respectively. 
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With these assignments to the two observed Pc(4380) and 
Pe(4450), their Jjil/p decay modes can be naturally inter¬ 
preted as the SfD* in (/ = 1/2,7 = 3/2) and the S*D* in 
{I - 112, J = 5/2); the molecular states can transit into Jlipp 
via exchanging an S -wave charmed meson. 

Let us investigate the decays further, where the Y,cD* state 
with (/ = 112, J - 3/2) and the S*Z)* state with (I - 112, J - 
5/2) transit into 7/(Ap. In the first process, the products 7/(A 
and p reside in an S-wave, whereas, for the second mode, be¬ 
cause the spin of 7’c(4450) is 5/2, the final 7/i/f and p must 
be in a D-wave to guarantee conservation of total angular mo¬ 
mentum. 

Usually, a D-wave decay is suppressed compared with an 
5-wave decay. Thus we can qualitatively explain why the 
width of f’c(4450) is much narrower than that of Pc(4380) 
[16]*. Additionally, a similar decay mode of these two hidden- 
charm molecular pentaquarks is rjcN, which is a D-wave de¬ 
cay mode, where N denotes a nucleon. 


a binding energy (E — -80 MeV) of the 'LcD* system with 
(/ = 3/2,7 = 1 /2) is the same as that of the ScD* system with 
(I - 1/2,7 = 3/2). In addition, we also find a bound-state 
solution for the system with (/ = 3/2,7 = 1/2) with 
a binding energy of -28 MeV if taking the same parameters 
as for the case of the system with (/= 1/2,7 = 5/2). 
Thus, there may exist two extra hidden-charm molecular pen¬ 
taquarks, the ScD* state with (/ = 3/2,7 = 1 /2) and the S*Z)* 
state with (/ = 3/2,7 = 1/2), which are the isospin parters of 
Pe(4380) and 7’c(4450), respectively . 

The experimental search for these two predicted isospin 
parters of ^^(4380) and ^^(4450) is an intriguing issue, 
which can be taken as a crucial test of the molecular as¬ 
signment of Pc(4380) and Pc(4450). Since the predicted two 
hidden-charm molecular pentaquarks are isospin-3/2 states, 
A(1232)7/i/^ and A(1232)77c can naturally be their decay prod¬ 
ucts. 


JcD* (1=1/2, J=3/2) 

501 - 



ZcD* (1=3/2, J=l/2) 



eJd* (1=1/2, J=5/2) 



EcD* (1=3/2, J=l/2) 



FIG. 1: The variations of the obtained OPE effective potentials for 
the systems to r, and obtained bound state solutions. Here, 

the masses of Pc(4380) and Pc(4450) can be reproduced well under 
the with (/ = 1/2, 7 = 3/2) and EJD* with (/ = 1/2, 7 = 5/2) 
molecular assignments, respectively. A = 2.35 GeV and A = 1.77 
GeV are taken for the and Z/5* systems, respectively. The 
blue curves are the effective potentials, and the red line stands for the 
corresponding energy levers. Additionally, the obtained spatial wave 
functions are given here. 

Besides explaining the observed Pc(4380) and Pc(4450), 
we further predict two hidden-charm molecular pentaquarks 
(see Fig. 1 (c) and (d)). We notice that the the OPE effective 
potential is the same for the EfZ)* system with {I - 1 /2,7 = 
3/2) and the l,cD* system with (/ = 2112, J - 1/2), and the 
only difference comes from their isospin and spin combina¬ 
tions. If taking the same parameters as input, we find that 


TABLE II: The typical values of the obtained bound state solutions 
[£(MeV), A(GeV]) for hidden-bottom and fi,-like zS.**B* and 
systems. 


(L7) 


^bD- 

ZfcB* 

(1/2,1/2) 

X 

X 


X 

(l/2,3/2) 

[-0.27, 1.22] 

[-0.26, 

1.34] 

[-0.27, 0.84] 


[-2.58, 1.32] 

[-2.62, 

1.44] 

[-2.36, 0.94] 


[-7.48, 1.42] 

[-7.63, 

1.54] 

[-6.88, 1.04] 

(3/2,1/2) 

[-0.27, 1.22] 

[-0.26, 

1.34] 

[-0.27, 0.84] 


[-2.58, 1.32] 

[-2.62, 

1.44] 

[-2.36, 0.94] 


[-7.48, 1.42] 

[-7.63, 

1.54] 

[-6.88, 1.04] 

(3/2,3/2) 

X 

X 


X 

(L7) 



E-fi- 

(1/2,1/2) 

X 

X 


X 

(1/2,3/2) 

X 

X 


X 

(l/2,5/2) 

[-0.28, 0.88] 

[-0.14, 

0.96] 

[-0.30, 0.64] 


[-3.18, 0.98] 

[-2.78, 

1.06] 

[-3.11,0.74] 


[-9.67, 1.08] 

[-8.97, 

1.16] 

[-9.51,0.84] 

(3/2,1/2) 

[-0.42, 1.02] 

[-0.30, 

1.12] 

[-0.28, 0.72] 


[-3.33, 1.12] 

[-3.03, 

1.22] 

[-2.74, 0.82] 


[-9.37, 1.22] 

[-8.91, 

1.32] 

[-8.19, 0.92] 

(3/2,3/2) 

X 

X 


[-0.28, 1.44] 


X 

X 


[-3.28, 1.60] 


X 

X 


[-9.13, 1.74] 

(3/2,5/2) 

X 

X 


X 


* The Jjif/p invariant mass spectmm around 4450 MeV is very complicated 
[16]. Another possibility is that Pc(4450) could be the P-wave excitation 
of the 5-wave molecular states which decays into JliJ/p mainly via 
P-wave. Hence its decay width is not very lai'ge. 


If the hidden-charm molecular pentaquarks, indeed, exist, 
there should also exist the hidden-bottom pentaquarks in ana¬ 
log to Pc(4380) and £^(4450). Based on the obtained OPE 
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effective potentials shown in Eqs. (8) and (9), we extend the 
same formalism to the B* pentaquark system. 

We need to specify that the hidden-charm and 

hidden-bottom have the same quantum numbers. Thus, 
the OPE effective potentials are completely the same. The 
reduced masses of the hidden-bottom molecular pentaquarks 
are larger than that of hidden-charm molecular pentaquarks, 
which means that the binding of B* should be more stable 
than that of the hidden charm pentaquark^. 

In the above discussion, we mainly focus on the exotic pen¬ 
taquarks which possess hidden charm or bottom. We can also 
extend the whole scenario to discuss the exotic states with 
open charm and bottom. The and systems which 

are Bc-like molecular pentaquarks may also exist. In our later 
works, we will carry out more research on the exotic states 
with open charm and bottom and make predictions on the 
spectra of the Bc-like and pentaquarks and their 

decay behaviors as well as the corresponding Bc-like mesons. 

According to the results presented in Table II, we draw our 
conclusions; 

1. There exist the S*B*, ScB*, bound states with ei¬ 
ther (/= 1/2, 7 = 3/2) or (7 = 3/2, 7= 1/2). TheS^B*, 
ZfB*, and 'LbD* states with (7 = 1/2,7 = 3/2) mainly 
decay into T(15')77/T(25')77, Bc(l^)N, and Bc(l“)77, re¬ 
spectively, while the typical decay modes of the S*B*, 
ZcB*, and S*7)* states with (7 = 3/2,7 = 1/2) include 
T(15)A(1232), Bc(r)A(1232), and Bc(r)A(1232), re¬ 
spectively. 

2. We can also find bound state solutions for the 2^B*, 

S*B*, S-wave systems with quantum numbers 

(7 = 1/2,7 = 3/2) and (7 = 3/2,7 = 1/2). In 
addition, the Z^B* state with (7 = 3/2,7 = 3/2) 
also exists. The main decay modes of S^B*, S*B*, 

with (7 = 1/2,7 = 3/2) are T(15)A/T(2B)A, 
Bc(l^)N, and Bc{l~)N, respectively. The S^B*, S*B*, 
2^77* states with (7 = 3/2,7 = 1/2) mainly decay into 
T( 15)A( 1232), Bc( 1 -)A( 1232), and Bd 1 -)A( 1232), re¬ 
spectively. T(1B)A(1232) is the main decay channel of 
the 2^B* state with (I -312, J - 3/2). 

In summary, the two newly observed resonant structures 
Bc(4380) and 7’e(4450) in the Jjil/p invariant mass spectrum 
of Kb JI xjjpK [16] are first identified as the hidden-charm 
molecular pentaquarks 2c7)* with (7 = 1/2,7 = 3/2) and 
2*77* with (7 = 1/2,7 = 5/2), respectively. Their mass spec¬ 
trum and qualitative decay behaviors are consistent with the 
existing experimental findings. 


^ Equivalently, the cutoff A in the form factor for the systems 

could be smaller than that for the systems if a solution for the bound 

states containing bb and cc of the same quantum numbers is reached. Ac¬ 
cording to experience for studying the S-wave systems, we try to 

search for the bound state solution in the range of A < 2.35 GeV for S^B* 
and the range of A < 1.77 GeV for 


The observation of Bc(4380) and 7’c(4450) has opened a 
new portal for investigating fermionic exotic states, i.e the 
long-searched mysterious pentaquarks. We propose that they 
are loosely bound molecular states composed of 2^*^ and 
Our study indicates that there should exist a 2^77* state with 
(7 = 3/2,7 = 1/2) and a 2*5* state with (7 = 3/2,7 = 1/2), 
which can be searched for via the A(1232)7/(A final state. The 
isospin partners of 7’c(4380) and 7’c(4450) shall be a crucial 
test of the present proposal. 

Under the present molecular assignments, the parity quan¬ 
tum numbers of both 7’£.(4380) and 7*^(4450) should be nega¬ 
tive. Although the present data of LHCb favor that 7’c(4380) 
and 7*^(4450) have opposite parities, they also mention in their 
paper that the same parities are not excluded [16]. It is highly 
probable that both states may be identified as possessing the 
same negative parity when more data and analysis are avail¬ 
able in the near future. If it really turns out that 7’c(4380) 
and 7’c(4450) have opposite parities, the physics will be very 
interesting. There are several options: (a) both of them are 
pentaquark states instead of loosely bound molecular baryons, 
as discussed in Ref. [30], (b) the lower state is an B-wave 
molecule while the higher state is a pentaquark state, or (c) 
the lower one is the S -wave molecule while the higher state 
is the T’-wave excitation. Then, one has to determine why the 
energy of the T’-wave excitation is as large as 70 MeV. 

Besides the above predictions, we also extend our formal- 
ism to the hidden-bottom 2)^ B* pentaquark and the Bc-like 
2/ B* and 2^ D* molecular pentaquark systems which con¬ 
tain open charm and bottom quarks. Several pentaquarks with 
hidden-bottom y}^^B* and the Bc-like pentaquarks and 

2^*^D*, which can be considered as the partners of 7’c(4380) 
and 7 ’c( 4450), are predicted. Experimental exploration of 
these predicted exotic states is a potential and important re¬ 
search topic for future LHCb experiments. 

Certainly, it is not the end of the story for the charmed and 
bottom pentaquarks. The discovery and theoretical studies on 
7’e(4380) and 7’c(4450) are just opening a new page of hadron 
physics. In the coming years, joint theoretical and experi¬ 
mental efforts will be helpful in pushing the relevant study 
on pentaquarks with heavy flavors, which will deepen our un¬ 
derstanding of non-perturbative QCD behavior further. 
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